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Abstract. We investigate the long term behavior for a class of competition-diffusion systems 
of Lotka-Volterra type for two competing species in the case of low regularity assumptions on 
the data. Due to the coupling that we consider the system cannot be reduced to a single equa- 
tion yielding uniform estimates with respect to the inter-specific competition rate parameter. 
Moreover, in the particular but meaningful case of initial data with disjoint support and Dirich- 
let boundary data which are time-independent, we prove that as the competition rate goes to 
infinity the solution converges, along with suitable sequences, to a spatially segregated state 
satisfying some variational inequalities. 
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<^ , 1. Introduction 

' Let Q be a bounded, open, connected subset of M. N with smooth boundary and let k be 

a positive parameter. The aim of this paper is to investigate the asymptotic behavior of a 
competition-diffusion system of Lotka-Volterra type for two competing species of population of 
densities u and v, with Dirichlet boundary conditions, 

Ut — Au = /(it) — kuv 2 , in f2 x (0, oo), 

Vt — Av = g(v) — kvu 2 , in x (0, oo), 

u(x, t) = ip(x, i), on dCl x [0, oo), 

v(x,t) = ((x,t), on dfi x [0, oo), 

u(x,0) = uq(x), in Q, 

v(x, 0) = vo(x), in 0. 

A relevant problem in population ecology is the understanding of the interactions between dif- 
ferent species, in particular in the case when the interactions are large and of competitive type. 
As the inter-specific parameter k ruling the mutual interaction of the species gets large, com- 
petitive reaction-diffusion systems are expected to approach a limiting configuration where the 
populations survive but exhibit disjoint habitats (cf. [7,8,10,19,22,24]). For population dynam- 
ics models which require Dirichlet boundary conditions we refer to [8,22], while for the more 
ecologically natural Neumann boundary conditions we refer to [10, 15] and references therein. 
As pointed out in [8], the Dirichlet case presents further difficulties compared with the Neu- 
mann case, as the boundary terms which pop up after integration by parts cannot be estimated 
independently of k. The classical stationary Lotka-Volterra model for two populations 

- kuv, in O, 
kvu, in Q, 
on <9$7, 
on dQ, 



(1.1) 
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has been intensively studied with respect to the spatial segregation limit as k — > oo. If, for 
instance, V an d C belong to W 1,OD (dQ), then there exists a sequence of solution (u K ,v K ) to 
(1.1), bounded in W 1,co (fl), and a limiting positive state (u,v) with uv = 0, satisfying suitable 
variational inequalities and such that, up to a subsequence, u K — > u and w K — > v in 
with a precise rate of convergence (see [5]). Concerning the parabolic system associated with 
(1.1), in [8] Crooks, Dancer, Hilhorst, Mimura and Ninomiya proved (also in the case of possibly 
different diffusion coefficients) that, for any T > 0, there exists subsequences u Km and v Km of the 
solutions converging in L 2 (fi x (0, T)) to a bounded state with disjoint support and solving a 
limiting free boundary problem. Beside this convergence results on finite time intervals, in [7], in 
the case of equal diffusion coefficients and stationary boundary conditions, Crooks, Dancer and 
Hilhorst recently studied the long term segregation for large interactions, by reducing the system 
to a single equation whose solutions admit uniform estimates in k. Typically, stabilization is 
based upon a variational structure yielding an energy functional, bounded and decreasing along 
the trajectories (see e.g. [13,27]). Unfortunately, as far as we know, due to the coupling term 
—kuv, the parabolic system associated with (1.1) does not admit a natural Lyapunov functional 
and a direct analysis is therefore not possible. Now, system (P K ) can be regarded as a variant 
of the standard Lotka-Volterra model, with different inter-specific competition coupling terms. 
In addition, if one considers homogeneous boundary data, then (P K ) admits a natural non- 
increasing energy functional A K : [0, oo) — > R 

A K (t) = -\\Vu(t)\\t + -\\Vv(t)\\l- / f(a)da- / g(a)da + - / u\t)v\t). 

As we will see, a non-increasing energy functional can be constructed also for general boundary 
conditions (see the proof of Theorem 2.11). We shall tackle the problem with techniques from 
the theory of dissipative dynamical systems to show the convergence towards the solutions to 
the stationary system, formally written as 

-Am = f(u) — kuv 2 , in S7, 

—Av = g(v) — kvu? , in 17, 

U = V'ooj on 

v = Coo, on dQ. 

A question which naturally arises is whether the solutions stabilize towards a segregated state 
along some tj — ► oo and kj — ► oo, for instance in the natural case when the initial data have 
disjoint support and the boundary data are stationary in time (see problem (1.2) in the next 
section). Some numerical computations in a square domain in R 2 have been performed in [8, see 
Sections 1 and 4] for the Lotka-Volterra model under these assumptions on the initial and 
boundary conditions (see also [17], where an algorithm for parallel computing was implemented 
in order to efficiently track the interfaces). In [26] we arranged a complete set of numerical 
experiments both for (1.2) (i.e. system (P K ) with time-independent boundary data) and the 
corresponding model with the standard Lotka-Volterra coupling. Although on one hand working 
with (P K ) gives some advantages in the study of the long term dynamics for k fixed as it directly 
admits a Lyapunov functional, on the other hand the asymptotic analysis for the solutions of 
(S K ) is far more complicated than the study of (1.1) (subtracting the equations of (1.1) one 
reduces to the single equation Au = ku{u — <I>) where $ is an harmonic function, while this is 
not the case working with (S K )). For instance, the global boundedness in k of the solutions in 
H 1 will be derived from the corresponding boundedness for the solution flow of the parabolic 
system uniformly with respect to k. To show the boundedness directly on the elliptic systems 
seems out of reach. In addition, the blow up analysis based on Lipschitz rescalings performed 
in [5] does not seem to work. 
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Concerning some physical motivations to consider coupling terms between the equations which 
are different from the standard one uv, we refer the reader, e.g., to Section 3.3 of classical 
Murray's book [21] (looking at formula (3.14) at p. 87, our system corresponds to the choice 
F(N,P) = 1 - N - kP 2 and G(N,P) = 1 - P - kN 2 with respect to the book's notations). 
It is also useful to think about systems of two Schrodinger [1] or Gross-Pitaevskii [9] equations 
modelling particle interaction (and populations can also be thought as discrete collections of 
interacting particles), intensively investigated in recent time (nonlinear optics, Bose-Einstein 
binary condensates, etc.), which present all the coupling of (P K ), yielding a variational structure. 
We refer the reader to [20] for the case k < 0, with physical motivations e.g. from [11], and to [23] 
for the case where k > 0, with physical motivations e.g. from [4]. Both [20,23] deal with the 
semi-classical regime analysis. 

1.1. The main result. The main result of the paper concerns with the long-term behaviour in 
large-competition regime for the system with time-independent boundary data, that is 

u t — Au = f(u) — kuv 2 , in 0, x (0, oo), 

vt — Av = g(v) — kvu 2 , in x (0, oo), 

u(x, t) = ip(x), on x [0, oo), 

v(x,t) = ((x), on 90 x [0,oo), 

u(x, 0) = uo(x), in f2, 

v(x, 0) = vo(x), in fi. 

Concerning the functions /, g : R — > R, let: 

(1.3) f,geC\[0,^)), f( s )=g(s) = 0, for all s < 0, 
f(s) < 0, g(s) < 0, for all s > 1, 

and we set 

F(t) = [ f(a)da, G(t) = [ g(a)da. 
Jo Jo 
The initial and boundary data are required to satisfy: 

(1.4) u ,v G H 1 ^), < uq{x) < 1, < v (x) < 1, a.e. in O,, 

(1.5) V,Ce# 1/2 (30), V = ^, C = v \an, 

(1.6) < ij>{x) < 1, < C(x) < 1, on dfl. 

Under these assumptions, as well as those of Section 2, for all k > 0, system (P K ) admits a unique 
global solution u K ,v K G C°([0, oo), H 1 ^)) n C 1 ((0, oo), L 2 ($7)). For the local existence, we refer 
the reader to a paper by Hoshino-Yamada [16] (see e.g. Theorems 1 and 2, having in mind to 
choose 6* = a = 7 = ^in Theorem l(i) and 7 = in Theorem 2(h), with respect to the notations 
therein). For smoothing effects we also wish to refer to the classical book of Henry [14]. The 
global existence result can be deduced by the comparison principle for parabolic equations (see, 
for example, the book of Smoller [25]). For u t — Au = f(u) — kuv 2 , vt — Av = g(v) — kvu 2 with 
positive initial data, one can show < u(t) < U(t) and < v (t) < V(t), where U, V are solutions 
of Ut — AU = f(U), Vt — AV = g{V) with the same initial and boundary conditions. Since U and 
V exist globally in time due to assumptions (1.3), (1.4) and (2.2) (a priori uniform-in-time L°°- 
estimates for the solutions hold, see Lemma 2.3), one also recovers the global existence result 
(for the sake of completeness, we also mention Theorem 3 in Hoshino-Yamada [16] for small 
initial data and part (iv) of Proposition 7.3.2 in [18] for smooth initial data). In the following 
we set M = H 1 ^) x H 1 ^), endowed with the standard Dirichlet norm, and 

H = {(u, v) G EI : uv = a.e. in O}. 
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The following is the main result of the paper, regarding system (1.2). 

Theorem 1.1. Assume (1.3) -(1.6) and (uo,Vo) 6 Ho. Let (u K ,v K ) be the solution to system 
(1.2). Then there exist two diverging sequences (ft m ), (t m ) C R + and (^ooi^oo) £ Hfe smc/i i/iai 

(n Km (t m ),f Km (t m )) -»■ (uoo^oo) in the V x L p norm, for any p £ [2,oo), 

as m -> oo, where 

Uoo,^oo>0, -AUoo < /(Uoo), -AVoo < p(Uoo), U^an = ?/>, «oo|an=C 

Moreover, in the one- dimensional case, we have 

\\(u Krn {t m ),v Krn (t m )) - {u 00 ,v 00 )\\ L c X > xL oo -> 0, as m ^ oo. 

Hence, starting with segregated data, the system evolves towards a limiting segregated state 
satisfying suitable variational inequalities. As we have previously pointed out, in Sections 1,4 
of [8], the reader can find very nice pictures reproducing (for the classical model) these kind of 
separation phenomena. Notice that, due to the nonstandard coupling in system (S K ) the H 1 
convergence seems pretty hard to obtain either working directly on the system (which would 
require precise quantitative estimate of the rate of convergence of the solutions to and i^) 
or using indirect arguments such combining blow up analysis with Liouville theorems (which, 
however, would naturally require stronger regularity assumptions on the boundary conditions). 
In Section 2, we will obtain, for k fixed, the asymptotic behaviour of the system in the case of 
almost stationary boundary data. The author is not aware of any other result of this type in 
the literature (see also [3]). 

2. Long term behaviour for k fixed 

The goal of this section is the study of the long term behaviour of the parabolic system (P K ), 
for any k > fixed. We cover the general case of boundary data depending on time. Finally, 
in the particular case of segregated initial data and time independent boundary conditions, we 
will prove a stronger global boundedness result. 

2.1. Assumptions and main result. Concerning / and g we will assume condition (1.3). 
Moreover, the initial and boundary data are required to satisfy (1.4) and 

(2.1) ^C€C°([0,oo), J ff 1 / 2 (an)) ) ^(0)=«o|an, C(0) = «o|an, 

(2.2) < ip(x,t) < 1, < C(x,t) < 1, on dtt x [0,oo). 
We will assume that: 

(2.3) ^MJ-^oo and CM) - * Coo in ff 1/2 (9fi) as t -> oo, 

(2.4) Vt,Ct G ^^oo-H^idn)) nL^CoojF- 1 ^^)), M;t)Xt(;t) -► as t -► oo, 

(2.5) Vt(-,0) =Ct(-,0) = 0, 

(2.6) ^CweL^Coo;^- 1 / 2 ^)). 

Under the previous assumptions we have the following result. 

Theorem 2.1. Let (uq,vq) G HI and k > 0. Then for every diverging sequence (t^) C M + there 
exist a subsequence (tj) C M + and a solution (u K ,v K ) £ M to system (S K ) such that 

\\(u K (tj),v K (tj)) - (u k ,v k )\\m -> 0, as j -> oo. 

Moreover, the convergence holds in the L v x L p norm for any p G [2, oo). 

Strenghtening the assumptions we obtain the global boundedness uniformly in k. 
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Theorem 2.2. Assume that (uo,vq) € Ho and the boundary conditions are time-independent. 
Then, in addition to the conclusion of Theorem 2.1, we have 

supsup||(u K (t),u K (t))|| H < oo, 

i>0 k>0 

namely (u K ,v K ) is bounded in H (and in any L p x L p space), uniformly with respect to k. 

This second achievement will be of course an important step in order to prove the main result 
of the paper. 

2.2. Some Preliminary results. From a direct computation, we have positivity and a priori 
bounds for the solutions to (P K ), uniformly with respect to k. 

Lemma 2.3. £ = [0, 1] x [0, 1] is a globally positively invariant region for system (P K ), uniformly 
with respect to k, namely 

< u K (x,t) < 1, < v K (x,t) < 1, a.e. x E CI, t>0. 

Proof. Testing the first equation of (P K ) with — u~ and using (1.3), (1.4), (2.1) and (2.2), we 
easily obtain that u K > 0, while testing the same equation with (u K — 1) + we deduce similarly 
that u K < 1. An analogous manipulation of the second equation in (P K ) yields the corresponding 
bounds for the component v K . □ 

Let A = —A be the Laplace operator on L 2 (VL) with domain S>{A) = Hq(Q) n H 2 (Cl) and 
consider the hierarchy of Hilbert spaces H a = 3)(A a l 2 ), a G R, with ||ii||_f/" = 1 1 ^4 CK / 2 ti 1 1 2 . We 
recall an exponential decay property of the heat kernel operator e* A . 

Lemma 2.4. Let a > 0. Then there exist uj > and C a > such that 

(2.7) \\e tA \\ C (L^H^) < C a e-"H- a , t > 0. 

In particular 



oo 

<tA| 



e CT \\c(L^,H^)d(T < oo, 





provided that a £ (0, 1). 



Proof. As the real part of the spectrum of A is bounded away from zero by a positive constant 
uj, by [14, Theorem 1.4.3, p. 26], for a > there exists C a > such that ||^4 a e" tj4 ||£( i 2 i 2) < 
C a e~ ult t- a , for all t > 0. Hence, \\e tA \\ c{L 2 tH 2 a) = ||(-A) a e* A || £(L 2 iL 2 ) < C Q e^~ a , for all 
t > 0. The second assertion follows by (2.7). □ 

Next we provide a compactness result for the trajectories of (P K ). 

Lemma 2.5. For any (u ,vo) £ M, k > and r > the set {(u K (t),v K (t)) : t > r} is relatively 
compact in M. 

Proof. Let U and V denote the solutions to the linear problems 

'U t -AU = 0, inftx(0,oo), 
U(x, t) = ip(x, t), on 90 x (0, oo), 
U{x,0) = U (x), inQ, 



(2.8) 
and 
(2.9) 



'V t -AV = 0, infix(0,oo), 
V(x,t) = £(x,t), on3Ox(0,oo), 
y(x,0) = V (x), infl, 
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where Uq, Vq € H 1 ^) satisfy 
i-AUo = 0, 



-AV = 0, in 9,, 

v (x) = C(x,o), onsn. 



in S7, 

I f7o(s) = -0(x,O), on Oft, 

By assumption (2.2) and the maximum principle for harmonic functions, < Uq(x) < 1 and < 
Vq(x) < 1 for a.e. Hence, arguing as in the proof of Lemma 2.3, we have < U(x, t) < 1 

and < V(x,t) < 1 for a.e. xe!l and i > 0. Now, the functions 

(2.10) u K (x, t) = u K (x, t) -U (x, t), v K (x, t) = v K (x, t) - V(x, t) 

solve the system with homogeneous boundary conditions 

' (u K ) t - Au K = f(u K + U)- k(u k + U)(v K + V) 2 , 

(v K ) t - Av K = g(v K + V)- k(v k + V)(u K + U) 2 , 

u K (x,t) = v K (x,t) = 0, 

u K (x,0) = uo(x) - U (x), 

v K (x,0) = v (x) - Vb(x), 

Denote now by * = ^>(x; t) G C°([0, oo), H 1 ^)) the family of harmonic extensions to O of tp 

j-A^(x;t) = 0, in Q, 
\ ^(x;t) = 4>(x,t), on dn, 



in x (0, oo), 
in f2 x (0, oo), 
on dQ, x [0, oo), 
in Q, 
in O,. 



(2.11) 



and set U(x,t) = U(x,t) — ^(x;t). Then U solves the nonautonomous problem with homoge- 
neous boundary and initial conditions 

'U t -AU = -^t, inftx(0,oo), 

(2.12) <U(x,t) = 0, on dn x (0,oo), 

U(x,0) = 0, inn. 

Notice that U(x,0) = since Uq(x) and $>(x;0) are both harmonic functions with the same 
boundary conditions. From (2.5)-(2.6) and classical regularity theory for harmonic functions, 

ll^t||L°°(0,oo;L 2 (n)) < c IIV't|lLoc»(0,oo;if- 1 /2(an)) ^ HV'tt Hz, 1 (0,oo;H^ 1 / 2 (dQ,)) 1 

By Duhamel's formula U is given by 

U(t) = - I' e (t ' a)A <f t {a)da. 
Jo 

If a e (1/2, 1), in light of (2.7) of Lemma 2.4, since ^ is in L°°(0, oo; L 2 (0)), 

(2.13) sup\\U(t)\\ H 2 a < oo. 

Of course the same control holds for V(t). Let now denote the harmonic extension of ipoo, 
the limit of ip(t) in H 1 ^ 2 (dfl) as t — ► oo according to (2.3) . By standard regularity estimates, 
\\^(t) — ^ooIIh^q) < c llV'(0 — V'ooIIh 1 / 2 ^)' so * na ^ ~^ in H 1 ^) as t — > oo. Of course 
the same control holds for the boundary extensions of (. Also, by Duhamel's formula we have 

u K (t) = e tA (u - U ) + f e^ A 1>l(a)da, 

Jo 

v K (t) = e tA (v - V Q ) + f e^ A <S> 2 K (a)da, 

Jo 

where 

$ «(o") = fMa)) - Ku K (a)v 2 K (a), $ 2 K (a) = g(v K (a)) - Kv K (a)u 2 K (a). 
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By means of Lemma 2.3, we have <I> 2 G L°°(0, oo; L°°(0)). If a € (1/2,1), then again by 
(2.7) one obtains for any r > 

(2.14) sup\\u K (t)\\ H 2 a < oo. 

t>T 

As H 2a is compactly embedded in H 1 ^) and u K (t) = u K (t) + U(t) + *(t) the assertion follows 
by (2.13)-(2.14) for the component u K . The same arguments works for v K . □ 

Remark 2.6. By strengthening the regularity assumptions on the boundary data, say W l, °°{dQ) 
in place of ^^{dVi) in the assumptions at the beginning of the section, and defining — A over 
L q {Vt) for any q > 2, the previous result can of course be improved, yielding compactness of the 
trajectories in W 2a,q (Q) for any q > 2, and hence into spaces of Holder continuous functions. 
Unfortunately the estimates are not independent of k and in order to have H 1 bounds uniformly 
in k we shall need to exploit energy arguments. 

For every r > and every function h : (0, oo) — > H 1 ^), let us set 

h T (t) = h(t + r), t>0. 

The following result gives a stabilization property for the solutions of the linear parabolic equa- 
tion with nonhomogeneous time-dependent boundary conditions. 

Lemma 2.7. Let U be the solution to the problem (2.8). Then U(t) — ► in ff 1 (f2) as t — ► oo, 
where Uqo G H 1 ^) is the solution to 



(2.15) 



—Alloc = 0, in fi, 
Uoo = ipoo, on dQ. 



Proof. With the notations introduced in the proof of Lemma 2.5, we consider, for r > 0, the 
functions W(t) = U T (t) - U(t) and g(t) = tf t (t) - *[(£), which satisfy 



lf(-A^ = ^), inftx(0,oo), 

W(x,t) = 0, on<9ftx(0,oo) 

W(x,0) = U{t) - C/ + *(0) - *(t), inn. 
By multiplying the equation by — AVF, we get 

d I|2 || A Tir/.\||2 



(2.16) 



^-IIVWWIIS + \\*w(t)\\i = - y n e(*)AW(t). 

By applying Holder and then Young inequalities on the right-hand side, we have 

j t \\vw(t)\\ 2 + \\\Awmt<\ut)\\l 

Let A be the positive operator on L 2 (f2) defined by A = —A, with domain T>(A) = H 2 (£i) n 
i?o(^)- Due to the (compact and dense) injection # 2 (ft)ni^(ft) = V(A) ^ V(A 1 / 2 ) = H^(Q), 

— 1/2 

we have ||VW||2 < ot\ ||AW||2 for some a± > (see e.g. Henry [14]), so that 

^\\vw(t)\\l + Y \\vw(t)\\ 2 <±\\ Q (t)\\ 2 . 

Finally, Gronwall inequality entails 

\\W(t)\\m < \\W(0)\\me- at + ce-^ f e° s \\ 6 (s)\\ 2 ds, 
Jo 

for some a > and c > 0. In turn, we readily obtain 

lim \\U T (t)-U(t)\\ m <4= lim ll*t(*) " *t(*)l|2 + Hm ||* T (t) - *(t)||ffi. 

t^oo \/0~ t— >oo t— »oo 
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In view of (2.3) and standard elliptic equations, we deduce 
(2.17) \\U T (t) - U{t)\\ H i -> 0, asi^oo. 

The same argument shows that {U(t)} t >o is bounded in H 1 ^). Let now (th) C M + be any 
diverging sequence. Since {U(t)} t >i is relatively compact in H 1 ^), there exists a subsequence, 
that we still denote by (th), such that U(th) — ► t4o in if 1 ^). Let rj G C^°(r2). By integrating 
the equation for £7 on + 1) xfl, yields 



lim / Utv+ / VU-Vrj = 0. 



On one hand, we have 
rt/,+1 



lim 

ft 



< lim / |tf(t h + 1) - tf(t h )|M < dim H^fa) - U(t h )\\ 2 = 0. 
ft Jn h 



it h Jn 

Moreover, there exists (sh) C K + with Sh = th + £h, < ^ < 1, such that by (2.17) 
th+i 



X7U-Vr]= / VU(s h )-Vi]= / VJ7(t fc ) • Vr? + o(l), as /i -► oo. 
t h ./n is Jn 

Hence, taking the limit as h — ► oo, we get J" n V£/oo • Vr/ = 0. Moreover, from the convergence of 
U(t h ) to £/oo in H 1 ^) we deduce that U(t h )\an -> C/oo|an in F 1/2 (<9$7). From (2.3) we deduce 
that [Too = t/>oo on <9$7. Therefore Uoo solves (2.15). Since (2.15) has a unique solution, we 
actually deduce the convergence of the whole flow U(t). □ 

Next, we obtain a summability result for the solutions to (2.8). 
Lemma 2.8. Let U be the solution to (2.8). Then U t G L 1 (0, oo; H 1 ^)). 

Proof. As in the proof of Lemma 2.5, U is the solution to (2.12). Hence, taking into account 
(2.5), it turns out that U(x,t) = Ut(x,t) is a solution to 



(2.18) 



'U t -AU = -V tt , inftx(0,oo), 
U(x,t)=0, on dn x (0,oo), 

U(x,0)=0, mil. 

By assumption (2.6) it follows * ti G L 1 (Q, oo; L 2 (Q)). In addition, we have * t G L 1 (0, oo; H 1 ^)). 
By Lemma 2.4 we have ||e* A ||_ C( - I/ 2 ( -^) ,h^(q)) — Ce~ wt t~ x l 2 , for some C, uj > 0. Hence, 

tf(t) = - fe^-^tti^da, 
Jo 

and we obtain 

rt 



l t/ 1lLi(0,oo;fl3(n)) ^ C 



-w(t— a) 



(t-o-r^WVttWhda 



o L Jo 

oo 



/"OO /'OO 

C/ ||*tt(a)|| 2 / e-^-^it-ay^dt 

JO l J a 

e -u>a a -l/2 d \u^ 



dt 
da 



C 



' tt\\ ^(0,00;^ (n))- 



Hence U G L 1 (0, 00; Hq(Q)), which yields Ut G ^ 1 (0, 00; Hq(Q)) and, in turn, taking into account 
(2.4), also U t G L^O, 00; H 1 ^)), concluding the proof. □ 
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Lemma 2.9. Let u K and v K be as in system P K . Then 

/ \\d t u K (a)\\lda < oo, / \\d t v K (a)\\l da < oo, 
Jo Jo 

for any T > 0. 

Proof. Setting T(x,t) = f(u K (x,t) + U(x,t)) - K,(u K (x,t) + U(x, t))(v K (x, t) + V(x,t)) 2 for any 
if!! and t > and = uo(^) — ^o(^)> it follows that u K is the solution to 

' dfU K — Au K = T in 0, x (0, oo), 

u K (x, t) = 0, on dd x [0, oo), 
^u K (x, 0) = m(x), in 17. 

Hence, since m G and T G L 2 (0, T, L 2 (17)) for any T > (as < u K ,v K < 1 and / is 

continuous), the desired summability for dtu K follows, e.g., by [12, Theorem 5, p. 360]. The proof 
for dtV K is similar. □ 

Let us recall a useful elementary Gronwall type inequality. 

Lemma 2.10. Let g G L^fO, oo), [0, oo)). Assume that T : [0, oo) — ► [0, oo) is an absolutely 
continuous function such that 

T(t) < ci + c 2 / g(a)^T(a)da, t > 0, 
Jo 

for some c±, c 2 > 0. Then 

T(t) < 2c 1 + cl||^||| 1(0;Oo) , t>0. 
Proof. Let t > and consider i G [0, i] such that T(t) = max{T(er) : a G [0, t]}. Hence 

T(t) < ci + c 2 / 5(cr)v / TMdcr < ci + c 2V / T(t)||5'||L 1 (o,i) < c i + C2V /T (^ll5 f llL 1 (0,oo)> 
Jo 

so the assertion immediately follows by Young inequality and Y(i) < T(i). □ 

Next we obtain an H 1 stabilization result for the solutions (u K ,v K ) to (P K ). 
Theorem 2.11. Assume that (uo,vq) G EI and set 

(2-19) = H^O^oIll + ll^i ll J L 1 (0,oo; J L 2 (^)) + || Z t \\ L i ( ,oo;L 2 (n)) • 

Then there exists a positive constant R = R(uo,vo,ip, C) independent of k such that 

(2.20) \\(u K (t),v K (t))\\n<R + Kn, forallt>0. 
Moreover, for any to > and k > 0, 

lim sup \\u K (t + t) - u K (t)\\ H i = 0, lim sup \\v K (t + t) - v K (t)\\ H i = 0. 

^°°re[0,r () ] ^°°Te[0,T„] 

Proof. Let tq > and k > 0. Let us first prove that 

(2.21) lim sup \\u K (t + t) - u K (t)\\ 2 = 0, lim sup \\v K (t + r) - v K (t)\\ 2 = 0. 

*^°°re[0,To] *^°°Te[0,To] 

According to the proof of Lemma 2.5, let again U (resp. V) be the solution of the linear problems 
(2.8) (resp. (2.9)), where Uq (resp. Vo) is the harmonic extensions of ^(0) (resp. C(0))- Then 
u K (x,t) = u K (x,t) — U(x,t) and v K (x,t) = v K (x,t) — V(x,t) are solutions to system (P K ) having 
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homogeneous boundary conditions. Let now e € (0, 1) and, taking into account Lemma 2.9, 
introduce the auxiliary energy functional A K : [0, oo) — > R defined by setting: 



K(t) = l\\Vu K (t)\\i + ^l|V5«(t)||l - J^F{u K {t) + Uit)) 
- ! G(v K (t) + V(t)) + ^ [ (u R (t) + U(t)f{v K {t) + V(t)f 

JQ 1 JQ 



+ 2 



+ 2 



Vu K (a) ■ VU t (a) 

f \ I Vv K (a) • VV t (a) 
Jo Uq 



o L JQ 

t 



du K (a) 



+ £ / ||^ti«(<7)|||d<T + e / ||^ K ((7)||lda. 

J ■/ 

We prove that A K is nonincreasing and there exist two constants a K £ R and /3 K £ R (which 
we will write down explicitely) such that a K < A K (t) < (3 K , for all t > 0. By multiplying the 
first equation of (P K ) by and the second one by dtv K , using the fact that U and V solve 
problems (2.8)-(2.9), and adding the resulting identities, we reaches 

(2.22) j t A K (t) = -(1 - e)\\d t u K {t)\\ 2 - (1 - e)\\d t v K (t)\\ 2 2 < 0- 
In particular {t A K (t)} is a nonincreasing function. Hence, 

A K (t)<A K (0) = hv(u -U ))\\ 2 2 + hv(v -V )g- [ F(uo)- [ G(v ) 

2 2 JQ JQ 

- ! W ■ V(n -U )- f W ■ V(v - V ) + £ [ u 2 v 2 , 

JQ JQ * JQ 

for alii > 0, namely A K is bounded from above, uniformly in time and f5 K is of the form 

(2.23) (3 K = P + k\\u v \\1 P = P(u ,v ,^,0- 

Now, using the trace inequality, the first equation of (P K ), the L°°-boundedness of the solutions 
and the Young inequality, we find c > and c £ > such that 



J°*£l,rM*)) da 



2||V* t (<7)|| 2 d<7+ / ||Afi K ((7)||2||*t(<7)||2d<7 
t 



[ ||V««((7)||2||V* t (<7)||2d7+ f 
JO JO 

f \\Vu K (a)\\ 2 \\Vy t (<r)hda + f \\dtu R (a)\\ 2 \\^ t (a)\\ 2 da 
Jo Jo 

f\\^ t {a)\\ 2 da 
Jo 

< [ \\Vu K (a)\\ 2 \\W t (a)\\ 2 da + e / \\d t u K {a)\\ 2 2 da 

JO JO 

+ c £ [ \\^t{a)\\ 2 2 da + CK f \\V t (a)\\ 2 da, 

JO JO 



< 



+ CK 



where is the harmonic extension of ipt to (see formula (2.11)). Analogously, we reach 
>dv K (a) 



o 



h-A dv 



'■M")) da < [ \\Vv K {a)\\ 2 \\VZ t {a)\\ 2 da + e [ \\d t v K (a)\\ 2 da 

/ JO JO 

+ c £ \\Z t (a)\\lda + ere / \\Z t (a)\\ 2 da, 
Jo Jo 
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where, instead, Z t denotes the harmonic extension of Ct to S7, namely 

j-AZ t (x;t) = 0, in Q, 

{ Z t (x;t) = Ct(M), on ^. 

From the above estimates, the definition of A K , (1.3), Lemma 2.7, and assumptions (2.4) we 
obtain that 

\\Vu K (t)\\ 2 2 + \\Vv K (t)\\ 2 2 <d + C 2 f \\Vu K (a)\\ 2 [\\VU t (a)\\ 2 + ||Wi(<7)|| 2 ]d<7 

Jo 

+ C 3 [ \\Vv K (a)\\ 2 [\\VV t (a)\\ 2 + \\VZ t (a)\\ 2 ]da, 
Jo 

for some positive constant C\ = C±(k) independent of t, 

(2.24) C 1 (n)=Q + K f i, Q = Q{u ,v ,ip,(), 

for C 2 ,Cs independent of t and n, where \i has been defined in (2.19). Hence, by the Cauchy- 
Schwarz inequality 

\\Vu K (t)\\ 2 + \\Vv K (t)\\ 2 <d + C 4 J* yfci+^ix 

x [||Vl7 t (<7)||2 + ||V* t (<7)|| 2 + ||W t (<7)|| 2 + ||VZ t (<7)|| 2 ] da 

for all t > 0, for some positive constant C4 independent of t and k. From Lemma 2.10 it follows 
that, for all t > 0, 

HVtkWIl! + ||V^(t)|| < 2d + C 4 2 [||W t || L 

+ ||VV^||£l( 0)O o;L 2 (n)) + ll V ^llL 1 (0,oo;L 2 (n))] 2 ) 

which, by Lemma 2.8 and assumption (2.4), yields boundedness of (u K (t),v K (t)) and, conse- 
quently of the sequence (u K (t),v K (t) in M, with the estimate appearing in (2.20). 

In particular, from the H boundedness of (u K (t),v K (t)) we deduce that A K is bounded from 
below uniformly with respect to t, with a constant a K of the same form as the one appearing in 
inequality (2.20) (say, A K (t) > — M — Nk/j,, for some constants M, N > 0). To prove this it suf- 
fices to repeat the estimates that we have obtained above (see the inequalities following formula 
(2.23)) on the term which appear in the functional as time integrals, using the H 1 bound of u K 
and v K , uniform in time. Notice that the time integrals e J * ||^t€t K (cr) |||rfcr and e J" * ||(?t£ K (cr)||2c/<7 
which appear in the estimate of the boundary term are balanced by the corresponding term in 
the definition of A K . More precisely, we obtain 
rt 



f f Vu K (a 
Jo Jn 



VU t {a)da 
VU(t) ■ Vu K (t) 



< 2{R + K/x)||VC/ t || L i (0iOO;i2(Q)) < A + BKfi, 
<(R + Kfi) sup \\VU(t)\\ 2 <C + Dkh, 



t>i 



as well as 



Jo H ~A du ,nK '/ H h Jo H ~A du ^ tK, / H h 

+ £ \\d t u K (a)\\l da + e \\d t u K {a)\\ 2 2 da > 
Jo Jo 

- ( J R + K/x)||V* t || i i (0iO o,L2 (n)) -C £ ||*t|||2 ( o )00iL 2 (Q)) -oe||*t|| £ i( 0) oo,La(n)) 

- (-R + K//)||VZ t || L i( 0)OOj L2(n)) - Ce||^t|li2(o,tx5,L 2 (n)) - CK ll^t|lLi(0,oo,L2(n)) > 



- E - Fkh, 
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for some constants A, B, C, D, E, F > independent of re and t. Now, for all r € [0, To], 

IK(t)-««(*)ll2 + IKW-««(*)lli 

= / |u^)-u K (t)| 2 + / \vl(t) - v K (t)\ 2 

rt+T rt+T 

<t \\d t u K (a)\\l da + t \\d t v K (a)\\ 2 2 da 
Jt Jt 

r t+T / d \ 

= Thj t ( " -^ K{<T) ) da ~^ ^ {t) " K{t + T0) ] ' 

where we exploited Holder inequality, Fubini's Theorem and identity (2.22) (in the spirit of [6]). 
Hence, we obtain 

(2.25) \\<(t)-u K (t)\\ 2 + \\v T K (t)-v K (t)\\ 2 

< 2(\\ul(t) - u K (t)\\ 2 + \K{t) - v K (t)\\ 2 + \\U T {t) - U(t)\\l + \\V T (t) - V(t)\\ 2 ) 

< ^ [A K (t) - A K (t + T )] + 2\\U r (t) - U{t)\\ 2 2 + 2\\V T (t) - V(t)\\l 

Since A K is nonincreasing and bounded from below at fixed k, it follows that A K (t) admits a 
finite limit as t — > oo. Therefore, letting t — > oo in (2.25), and taking into account Lemma 2.7, 
we obtain (2.21). Now, assume by contradiction that, for some e$ > 0, 

\\u K (t h + T h ) - u K (t h )\\ H i > e > 0, 

along a diverging sequence (th) C M + and for (t^) C M + bounded. In light of Lemma 2.5, there 
exist u and u 6 i^ 1 such that, up to a subsequence that we still denote by (th), u K (th+Th) — > n 
in i? 1 (ri) as /i — > oo, and u K (th) —> u in // 1 (0) as /i — > oo. In particular, ||n — > £o > 0, 

while (2.21) yields \\u — u\\l2 = 0, thus giving rise to a contradiction. One argues similarly for 
v K . This concludes the proof of the theorem. □ 

Next we have an important consequence of the previous lemma, proving Theorem 2.2. 

Corollary 2.12. Assume that (uq,vq) £ Mq and that the boundary data are stationary. Then 
the sequence (u K (t),v K (t)) is uniformly bounded in H l with respect to t and k. Moreover the 
energy functional which appears in the proof of Theorem 2.11 is bounded below and above by 
constants which are independent of k. 

Proof. If (no, vq) £ Ho, since no^o = and ipt = Ct = by (2.19) we have that fi = 0. In turn, by 
(2.20), the sequence (u K (t),v K (t)) is uniformly bounded with respect to t and k. By inspecting 
the proof of Theorem 2.11 it is easy to check that the auxiliary energy functional satisfies 

—M - Nk^l < A K (t) <0 + Pkh, t > 0, 

for some constants M, N, O, P > independent of k. Hence, being \x = it follows that A K has 
bounds uniform in time and in k. □ 

2.3. Proof of Theorem 2.1 concluded. Let re > and let (th) C M + be any diverging 
sequence. Then, by virtue of Theorem 2.11, we have 

(2.26) lim \\u K (t h + r h ) - u K (t h )\\ H i = 0, lim \\v K (t h + r h ) - v K (t h )\\ H i = 0, 

h^oo h— too 

for every sequence (t/J C [0, 1]. Let us fix r], £ G C£°(Q,). By integrating over (th,th + 1) x fi 
the equations of (P K ) multiplied by r\ and £ respectively, we reach 

rt h +l r rt h +l p rt h +l r rt h +l 



lim 

h 



\ / d t u K ri + / Vu K ■ Vr? - / / f(u K )rj + re / / u K v 2 K n = 0, 

Jt h Jn Jt h Jn Jt h Jn Jt h Jn 
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lim 

h 



/ / d t v K £ + / / Vv K ■ V£ - / / g(v K )Z + k / v K u 2 K £ = 
Jt h Jn Jt h Jn Jt h Jn Jt h Jn 



Regarding the first terms in the previous identities, we obtain 



lim 

h 



lim 

h 



t h Jn 

t h +i 



d t u K rj 



n 



< lim / \u K (t h ) - u K (t h + l)\\<q\ < limc\\u K (t h ) - u K (t h + 1)|| 2 = 0, 

h Jn h 

< lim / \v K (t h ) - v K (t h + 1)| iCl < hmc\\v K (t h ) - v K (t h + 1)|| 2 = 0. 

h Jn h 



Moreover, there exist two sequences (sh), (?7i) C R + such that 

th< s h <t h + 1, t h <r h <t h + 1, s h = t h + p\, r h = t h + p 2 h , 
with (p\),(p 2 h ) C [0,1], and 

/ / Vu K -Vr]- f(u K )r] + nu K vlr] = / Vu K (s h ) ■ Vr] - f(u K (s h ))r] + KU K (s h )v 2 K (s h )f], 

Jt h Jn Jn 

! h ! [Vv K • V£ - g{v K )t + kv k u 2 £] = f Vv K (r h ) • V£ - g(v K (r h ))£ + KV K (r h )u 2 K (r h )£. 

Jt h Jn Jn 

In turn, we get 

/ Vu K (s h )-Vr]- / f(u K (s h ))r] + k u K (s h )v 2 K (s h )r] =0, 
Jn Jn Jn 

[ Vv K (r h )-V£- [ g(v K (r h ))£ + K [ v K (r h )u 2 K (r h )Z = 0. 
Jn Jn Jn 

On the other hand, in light of (2.26), there holds 

lim \\Vu K (s h ) - Vu K (t h )\\ 2 = lim \\Vu K (t h + p\) - Vu K {t h )\\ 2 = 0, 

h h 

\im\\Vv K (r h ) - Vv K (t h )\\ 2 = lim \\Vv K (t h + p\) - Vv K (t h )\\ 2 = 0. 

h h 

Hence /^^^(s/,) — Vu K (^)) • V?? — ► and, as f,g are C 1 on [0, 1] and < u K ,v K < 1, 



lim 

h 



lim 

h 



(f(u K (s h ))-f(u K (t h )))r) 
as h — > oo, and, finally, 



< csup|/ / |||n K (s fe ) - u K (t h )\\ 2 -> 0, 
[o,i] 



(««(sft)«K(sh) - u K {t h )v 2 K {t h ))r] 



< c\\u K (s h ) - u K (t h )\\ 2 + c\\v K (s h ) - v K (t h )\\ 2 -> 0, 



as /i — ► oo, the positive constant c varying from line to line. Of course, the same conclusions 
hold for the limit involving the sequence v K (rh)- In conclusion, we reach 



lim 

h 



0, 



/ Vu K (t h )-Vrj- / f{u K {t h ))rj + k I u K (t h )vl(t h )rj 
Jn Jn Jn 

lim I Vv K (t h )-V£- ( g{v K {t h ))i + K ! v K {t h )u 2 K {t h )i =0. 
h VJn Jn Jn 

Again in view of Theorem 2.11, we can assume that, up to a subsequence, which we shall denote 

again by th, we have that u K (th) — 1 u K and v K (th) — v K weakly in H l {Q). Up to a subsequence, 

in light of Lemma 2.5, this convergence is actually strong. Notice also that 

u K \ d n = limu K (t h )\dn = lirmp{t h )\dn = ipoo, 

h h 

v K \an = timv K (t h )\dn = lim^(t h )| an = Coo, 

h h 
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where we exploited the compact embedding H l / 2 (dVl). Moreover, by Lemma 2.3 and 

the Dominated Convergence Theorem, as h —>■ oo, we get 

f Vu K -Vn- [ f(u K ) V + K [ u K v 2 K r] = 0, \/rj G Hq(£1), 
Jn Jn Jn 

f Vv K • V£ - / g(v K )£ + k [ v K uU = 0, V£ G H^Sl). 
Jn Jn Jn 

Hence (u K ,v K ) G H is a solution to (<S K ). The convergence occurs of course in L p (f2) for any 

p G [2,2*). For p > 2*, taking e > and using the bounds < u K (t h ) < 1 and < u K < 1, we 

have 



/ 



|««(tft)-««r < 2 p+£ - 2 *||n K (^) -u«||£:*, 



concluding the proof. □ 



3. Proof of Theorem 1.1 

Before concluding the proof of Theorem 1.1, we provide the convergence of the sequences 
(u K ,v K ) in any L p space with p > 2 towards a segregated state. Notice that the solutions 
(u K ,v K ) to (S K ) pop up as H 1 limits of the solutions to (1.2), and the boundedness of (u K ,v K ) 
in H 1 is inherited by the boundedness of (u K (th),v K (th)) in H 1 uniform in t and k (in the 
case (uq,vo) G Ho). Without this information it would not have been clear how to show the 
boundedness of (u K ,v K ) working directly on the elliptic system (instead, for system (1.1), this is 
an easy task, cf. [5, Lemma 2.1]). 

Lemma 3.1. Assume that (uq,vq) G Ho- Let (u K ,v K ) G HI be the solution to system (S K ) as 
obtained in Theorem 2.1 for k > 0. Then there exists (iioo,^oo) G Ho with 

■Uoo,Woo>0, -AUoo < /(Uoo), -A-Uoo < ^(Voo) 

and Uoolgn = ip, v^lgn = C such that, up to a subsequence, as k — > oo, 

(u K , v K ) — > (uqo, foo) in the L p x L p norm for any p G [2, oo). 

Proof. By virtue of Corollary 2.12 the sequence (u K (th),v K (th)) is bounded in H, uniformly with 
respect to k. Hence, since (u K ,v K ) is the iJMimit of (u K (th),v K (th)) as h — > oo, we deduce that 
(u K ,v K ) is bounded in H and < u K (x) < 1, < w K (x) < 1, for a.e. x G fi. Taking into account 
that some terms in the functional A K introduced within the proof of Theorem 2.11 vanish under 
the current assumptions (stationary boundary conditions) and that the terms e J * ||<9i-u K (<7)||| 
and e J * ||<9 t i; K (<7)||| were artificially attached to make things work (notice that the original A K 
is decreasing also in the case e = 0, see formula (2.22)), we now just consider the natural energy 
functional (for the sake of simplicity we do not change the name) 

A«(t) = ^||V««(t)||l + \\\Vv K (t)\\ 2 2 ~ J F(u K (t) + U(t)) 

- [ G (v K {t) + V{t)) + ± f (u K (t) + U(t)f{v K (t) + V(t)f. 
Jn 1 Jn 

Then, we have 

« f ««(tfc)^(«h)=2A«(t h )-||V««(tfc)||l-||Vi5«(tft)||l + 2 f F(u K (t h )) + G(v K {t h )). 
Jn Jn 

Since by Corollary 2.12 the right hand side is uniformly bounded with respect to k, we have 

(3.1) k [ ulvl < C, 

Jn 
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for some positive constant C independent of k. Let G and Voo G H 1 ^) be the weak 

limits, as k — > oo, of u K and v K in H l (Vl) respectively. Or course, by the compact embedding 
up to a further subsequence, u K — > and £ K — > in L P (Q) for any 
p G [2, 2*) and < Moo(^) < 1, < Woo(^) < 1, for a.e. i£H. In the case p > 2*, let e > 0, so 
that 

/ ~~ ^oo| P ^ 2 P ^ — U 00 ||2*_ e , 

yielding again the convergence. Due to inequality (3.1), we get 

lim / u 2 K v 2 K = / u^v^ = 0, 

which yields Uoo^oo = a.e. in Q, namely («ooj^oo) G Ho- Moreover, for each k > 0, 

-Au K < /(tt«), -Aw K < 

which pass to the weak the limit, yielding — A-Uqo < /(itoo) and — Afoo < g(voo). By the compact 
embedding i7 1 (0) H l / 2 (dfl), also the boundary conditions are conserved. □ 

3.1. Proof of Theorem 1.1 concluded. We can now conclude the proof of Theorem 1.1. Let 
(uq,vq) G Ho, p G [2, oo) and let (i/J C M + be any diverging sequence. In light of Theorem 2.1, 
for every k > 1, there exist a solution (u K ,v K ) of (S K ) and a subsequence (t%) C IR + such that, 

\\{u K (t%),v K (t%)) - (u K ,v K )\\m -»■ 0, as /i -»■ oo. 

Moreover, by Lemma 3.1, there exists (t/.oo, Woo) G H with the required properties, such that, 
up to a subsequence, 

) - (woo,Woo)||lpxLp — > 0, as K ► oo. 
Now, let m > 1 and let K m > 1 be such that 

II ("Km. "Km ) ~ (^OO , ^OO ) 1 1 LP X U> < ^ • 

Then, there exists t K h m > 1 such that 

\\(u Km (t K h 2),V Km (t K h 2)) - (u Km ,V K J\\ L rxLr < 

In turn, setting t m = t 1 ^™, and combining the previous inequalities, we get 

\\(u Km (t m ),v Km (t m )) - (u 00 ,w 00 )|| LPxLP < 

which concludes the proof of the first assertion. In the one dimensional case, by means of Morrey 
Theorem, for every x,y G ft, we have 

\u Km {t m ){x) - u Km (t m )(y)\ < A\\Vu Km (t m )\\2y/\ x -y\ ^ Cy/\x-y\, 

together with \u Km (t m )(x)\ < 1, yielding the convergence to («oo,^oo) in the L°° x L°° norm via 
Ascoli's Theorem. □ 
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